Quantum Lubrication: Suppression of Friction in a First Principle 

Four Stroke Heat Engine. 

Tova Feldmann and Ronnie Kosloff 1 

1 Department of Physical Chemistry, 
Hebrew University, Jerusalem 91904, Israel 
(Dated: February 9, 2008) 

Abstract 

A quantum model of a heat engine resembling the Otto cycle is employed to explore strategies to 
suppress frictional losses. These losses are caused by the inability of the engine's working medium 
to follow adiabatically the change in the Hamiltonian during the expansion and compression stages. 
By adding external noise to the engine frictional losses can be suppressed. 
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I. INTRODUCTION 



Working conditions of real heat engines are far from the ideal reversible limit. Their per- 
formance is restricted by irreversible losses due to heat transport, heat leaks and also friction. 
Actual working devices tend to optimize the performance by balancing the losses with maxi- 
mizing work|l|. For engines producing finite power irreversible losses are unavoidable. High 
performance engines are therefore constructed from materials which reduce heat resistivity 
while minimizing heat leaks. In addition lubricants are employed to reduce frictional losses. 
The present study explores quantum lubricants, schemes to reduce the frictional irreversible 
losses and thus enhance the performance of the quantum heat engine. 

Quantum models of heat engines based on first principles are remarkably similar to their 
macroscopic counterparts . These engines extract heat from a hot bath of temperature 



T h and eject heat to a cold bath of temperature T c . The irreversib 



c 



osses due to the finite 



rate of heat transport have been linked to their quantum origin jj, ^, (], Q]. Optimal perfor- 
mance strategies lead to solutions where the working fluid never reaches thermodynamical 
equilibrium with the heat baths. Performance curves can be directly compared to those ob- 
tained in finite time thermodynamics which employ phenomenological heat transport laws 

Friction is the punishment for compressing or expanding the working medium too fast. 
In a quantum engine, compression/expansion is a change in an external field described by a 
parametrically time dependent Hamiltonian of the working medium. Whenever the control 
Hamiltonian does not commute with the internal Hamiltonian of the working medium, the 
rapid change in the external field does not allow the state of the working medium to follow 
adiabatically the instantaneous energy levels [jjlQ)- As a result both coherences and 
additional energy becomes stored in the working medium. The dissipation of this additional 
energy in the cold bath together with the inevitable decoherence is the quantum analogue 
of friction. The key to quantum lubrication is to suppress the creation of off diagonal terms 
in the energy representation. 

The quantum four stroke Otto cycle is chosen to demonstrate the lubrication effect. 
The working medium is composed from interacting two-level systems. Accordingly, the 
uncontrolled internal Hamiltonian becomes: 

H mf = 2-- i / 2 j(*l®*l-* l y ®d- 2 y ) = JB 2 . (1) 



2 



where <x re pres ent the spin Pauli operators and J scales the strength of the inter-particle 
interaction jlp[llll|. The external control Hamiltonian is chosen as: 

H ext = 2- z ' 2 uj{t) (&l ® f + I 1 ® afj = uj(t)B! , (2) 

where ui{t) represents the external field. The total Hamiltonian becomes: 

H = cj(t)Bi + JB 2 , (3) 



where Q(t) = \J lo 2 + J 2 defines the temporary energy scale. At various times H(t) does not 
commute with itself since [Bi,B 2 ] = >/2iB 3 ^ , ( B 3 = 2~ 3/2 (&l ® & 2 X + cr\ ® d-J)). 
The set of operators {B} forms a closed orthogonal Lie algebra. In addition, ^B^ ■ Bjj = 
trjB^Bj} = and tr{Bj,} = 0. The irreversible equations of motion for this set are 
= i[H, B] + £^(B) where is the dissipative Liouville superoperator and the set 
will be also closed to C* D . A thermodynamical description requires that the set of variables 
{b(£)} = {(B)} should be closed to the dynamics generated on all branches of the engine's 
cycle. 

The energy balance of the engine is composed of the heat flow and power: 

n= v+Q < 4 > 

where: V = {% = li(fii> , and Q = (C D (H)). 

The state of the working medium p can be reconstructed from five thermodynamical 
variables 6k = (B&) composed of the expectation of the three operators in the Lie algebra 
and two additional ones, B 4 = 2~ 3 / 2 (&l <g> | 2 — f 1 <g> a 2 ^j and B 5 = \cr\ ® a 2 z [ll] leading 
to: 

P = ^I + X>B* • (5) 

k 

The occupation probability p n of the energy level n defines the energy entropy: 



S E = -^p„lnp n . (6) 

n 

If [p, H] then this entropy is different from the von Neumann entropy 

S = -tr{p\np} , (7) 
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and Se > <S- The difference between Se and S is a signature of friction The external 
entropy production is a measure of the irreversible dissipation to the hot and cold baths: 



where <2/j and Q c are the heat dissipated to the hot or cold baths respectively. 



II. THE CYCLE OF OPERATION 

A four stroke cycle of operation is studied. As shown in Fig. QJthis cycle includes: 

• An adiabatic expansion branch where an external field is chosen to decrease linearly 
from u>b to uj a 

• A cold isochoric branch where heat is transferred from the working medium to the 
cold bath (T c ). 

• An adiabatic compression branch where an external field is increased linearly from uj a 

to UJf). 

• A hot isochoric branch where heat is transferred from the hot bath at temperature 
to the working medium. 

This cycle is a quantum model of the macroscopic Otto cycle. The control parameters are 
the time allocations on the different branches, the total cycle time and the extreme values 
of the external field. 

The cycle of the engine becomes a sequence of four completely positive maps that de- 
fine the different branches. Eventually this sequence closes upon itself. Repetition of the 



121. The 



sequence of controls leads to steady state operational conditions or a limit cycle 
map Uk relates the initial set of these operators to their final values for each of the engine 
branches. These maps are obtained by solving the equations of motion for the set of op- 
erators {B}. The overall cycle map is the product of the individual maps of each branch 
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On the isochores the maps Uh/ C ar e generated by the completely positive generator C* D 
jlHI ]. Forcing detailed balance C* D leads to thermal equilibrium with a rate determined 
by T. In addition C* D also degrades the off diagonal elements of p e , interpreted either 
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FIG. 1: Three cycles of the heat engine in the (uj,Se) plane. The common values are : J = 
2., T h = 7.5, T c = 1.5, uj a = 5.08364, to b = 12.6355, T h = T c = 1.16748, -f h = -0.05, <y c = 
—0.06. Afe a = A ao = for the black cycle, A ba = 1.28, A ab = 0.64 for the red cycle , and 
Ab a = 122.88, A ao = 61.44. for the green cycle. The time allocations on the different branches 
correspond to the optimal power point, when A& a = A a & = 0. The optimal time allocations are: 
r h = 1.0795 T ba = 0.01478, r c = 1.0088 , r ab = 0.0069. The points E and F represent the 
equilibrium points of the hot and cold isochores. 



as decoherence or as dephasing. The dephasing time T2 becomes identical to the energy 
equilibration time T 2 = T\. The dissipation has to also eliminate the additional energy 
accumulated on the adiabat. Degrading the coherences causes the frictional process to 



become irreversible 



11 



The interaction of the wo riving medium with the bath can 
also be elastic. These encounters will scramble the phases conjugate to the energy, and the 
associated decay time is termed pure dephasing (T 2 *). In Lindblad's formulation it becomes 
C* D e(B) = — 7[H, [H, B]] and T 2 * = l/2'~fQ 2 . Note that elastic medium can not transfer or 
absorb heat, therefore we always need an inelastic medium. 

On the adiabats, the varying field u(t) causes an explicit time dependence of H 
= z[H(t),B]. Since the energy eigenvalues change, even if initially [p, H] = 
the state p will develop off-diagonal terms in the energy frame (Cf. Eq. (B6) in Ref. 
The external power V = ubi can also be decomposed to the diagonal and off-diagonal terms 
in the energy representation: 



(9) 
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The first diagonal term represents the power required to compress or decompress the working 
fluid: 

r f J M = ^ , (10) 

the second term in Eq. is the additional power required to drive the working fluid in a 
finite rate: 

VlT Uon = 6 ^((jB 1 -coB 2 )) . (11) 

This term represents the power invested against friction therefore it vanishes when J = or 
u = 1141. 



III. QUANTUM LUBRICATION 

A good lubricant should be able to increase the overall optimal power of the engine. The 
insight that energy coherences leads to frictional losses, suggests that forcing the cycle tra- 
jectory to follow adiabatically the instantaneous energy levels will be beneficial. A possible 
approach is to increase the dephasing on the isochores so that at the beginning of the adiabat 
[p, H] = 0. A thermal bath with strong dephasing will cause such an effect. We have found 
that an addition of pure dephasing has only a minor effect on the performance of the engine. 

The quantum "lubricant" has to suppress the creation of the energy coherences on the 
adiabats. Formally this can be described by a generator of dephasing in the equations of 
motion for the set {B} on the adiabat 

^ = i[H, B] - A[H,[H,B]] (12) 

and A is the dephasing coefficient. 

The success of this approach is shown in Fig. El As a reference the optimal power of the 
engine as a function of cycle time is shown around the global maximum. Each point on the 
graph is optimized with respect to the time allocations on the four branches of the cycle. 
Employing these time allocations the power of the engine is recalculated with the addition 
of the dephasing term on the adiabats. It is clear in Fig. Elthat in the interval of cycle times 
around the maximum power the "lubricated" engine outperforms the optimal solutions of 
the reference engine. The "lubricated" maximum power point also moves to shorter cycle 
times. 
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FIG. 2: The optimal power as function of cycle time r with and without the quantum "lubricant". 
The optimal time allocations are obtained for A;, a = A a t = (blue line). The red line is the power 
with Ah a = l-28A a & = 0.64. All the other parameters are as in Fig. ^ 



For longer time allocations on the adiabats where less external power is consumed to 
overcome the friction, the performance enhancement due to dephasing decreased eventually 
leading to a crossover where dephasing on the adiabats decreased the power. For larger J 
values we also found that dephasing was not able to improve the performance. 

Fig. H2 shows the accumulated work against friction y\]f nctlon — J pfncuon^ _ 
f M^-^ dt (Cf. Eq. (HU) ) function of time on the adiabat for increasing dephasing 
parameter. The main point is that increasing dephasing eliminates the work against friction. 
This improvement saturates once yyf rictlon i s eliminated. 

Another consequence of the quantum lubrication is that the energy entropy Se does 
not increase on the adiabats as can be seen in Fig. |I] Asa result the energy entropy Se 
approaches the von Neumann entropy S, Eq. (JJJ). These results establish the principle 
of quantum lubrication maintaining the working fluid in a diagonal state in the energy 
representation. 



A. Dephasing synthesis 

Suppression of friction requires a method to synthesize dephasing on the adiabats. The 
dynamics of the working medium has to be changed from unitary to dissipative. The obvious 
approach of adding a dissipative bath on the adiabatic branches is difficult to achieve. Such a 
bath should have only elastic encounters with a system with a time dependent Hamiltonian. 
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FIG. 3: The work performed by the engine on the adiabat b — > a branch, separated into wf rtcUon 
(left scale) and W* (right scale) for different dephasing values as a function of n& a (i). W^ teld 
has only a weak dependence on the dephasing value. The cycle parameter values are as in Fig. ^ 

The solution is to employ the external controls of the engine to synthesize the dissipation. 
The idea comes from the singular bath limit, a bath generated from a system operator 
coupled to a delta correlated noise H sr = A s g(t) where (g(t)g(t')) r = j5(t — t') where the 
average is taken over the bath fluctuations. The Liouville generator associated with this 
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3- 



system bath coupling becomes: £*(X) = — ^- [A s , [A s , X] 

To implement such a scheme random noise is added to the external controls of the engine. 
The implementation divides the adiabat branch into N segments. In each of these segments, 



the external field u is constant and is chosen to be (for the a — > b adiabat): oj^ = tu a + 



N 



k 



for the kth segment. The short time propagator on the k segment for the set B becomes: 



Z4B 



e «H(w fc )At fc g e ~iU(uj k )At k 



(13) 



where At*, is the time interval of the kth segment. At this point random noise is added to 
the time interval 



At k 



Tab 

N 



1 + r) 



(14) 



where r is a random number with zero mean and variance a. Expanding the propagator Eq. 
fH3j) to second order and averaging over the random noise will lead to the average generator 



for the k time segment: C* k (B) = i[H(o;^),B] 



Na 2 



[H(co>fc), [H(a>fc), B]]. In the limit when 



N — > oo this average propagator becomes identical to Eq. (fT2|) provided a 



2r ad A 
N ' 



The addition of random noise means that the individual cycle has to be replaced by 
the average performance on many cycles. As a result only an average cycle time can be 



S 




FIG. 4: Power (right scale) and entropy production (left scale) as a function of the dephasing 

1/2 

parameter AJ b and the variance, a a b of the random fluctuations of the time segment on the adiabats 
(upper scale). The dephasing synthesis results are represented by stars and filled triangles. The 
empty circles and squares represent the dephasing dynamics. The calculations were performed 
with 4«fc = Sat. 

A ba °"6a 

defined. This noisy lubrication procedure was simulated with N = 200 on both adiabats. 
The power and other thermodynamic variables were calculated as an average of 2000 cycles. 
Convergence was checked by continuing this averaging 1000 additional times. 

Fig. 0]compares the Power V and entropy production per cycle, AS u /t calculated by the 
two methods for the time allocations of the maximum power point. It is clear that the results 
obtained by imposing dephasing on the adiabats Eq. (fT2"|) are identical to the dephasing 
synthesis Eq. (|13|). The signature of lubrication is the reduction of entropy production 
which accompanies the increase in power. This is contrary to optimizing the power with 
respect to heat transport. In that case the increase in output power is accompanied by an 
increase in entropy production Q|. The choice of the procedure to generate dephasing is 
unique. For example, adding the random noise to the frequency u>k at each time segment 
has been tested. The performance of the engine only became worse. The reason is that such 
a term leads to the dissipative generator £d(A) = — ^[Bi, [Bi, A]] which does not eliminate 
the off diagonal elements in the energy representation. 

The present study should be related to other recent works. For example, adding me- 
chanical noise to a quantum refrigerator has been shown experimentally to cool atoms in a 



magnetic trap 



17| . It seems that the mechanism involves inducing non unitary dynamics. 
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Contrary to the present study, in ot 
lating the second law Scully et. al. 



: rer scenarios coherence can be beneficial. Without vio- 
3] showed that additional work can be extracted from 
the coherences in quantum heat engine. 

To summarize, frictional losses are caused whenever \H. ext , H in4 ] 7^ 0. Then the fast 
dynamics induces coherences in the energy frame. The essence of quantum lubrication is 
suppressing the generation of these off diagonal elements in the energy representation. The 
present model demonstrates how externally induced noise can achieve this task. 

Acknowledgments 

Work supported by the Israel Science foundation. We want to thank Lajos Diosi David 
Tannor and Peter Salamon for many helpful discussions and Moshe Goldstein for the random 
number generator. 



[1] P. Salamon, J.D. Nulton, G. Siragusa, T.R. Andersen and A. Limon, Energy 26, 307 (2001). 

[2] T. D. Kieu, Phys.Rev.Lett. 93, 140403 (2004). 

[3] K. H. Hoffmann, Ann. Der. Phys. 10, 79 (2002). 

[4] Eitan Geva and Ronnie Kosloff, J. Chem. Phys. 96, 3054 (1992). 

[5] Eitan Geva and Ronnie Kosloff, J. Chem. Phys. 97, 4398 (1992). 

[6] Tova Feldmann, Eitan Geva, Ronnie Kosloff and Peter Salamon, Am. J. Phys. 64, 485 (1996). 

[7] Eitan Geva and Ronnie Kosloff, J. Chem. Phys. 104, 7681 (1996). 

[8] F. Curzon and B. Ahlborn, Am. J. Phys. 43, 22 (1975). 

[9] B. A. P. Salamon, A. Nitzan and R. S. Berry, Phys. Rev. A 21, 2115 (1980). 

[10] Ronnie Kosloff and Tova Feldmann, Phys. Rev. E 65, 055102 1 (2002). 

[11] Tova Feldmann and Ronnie Kosloff, Phys. Rev. E 68, 016101 (2003). 

[12] Tova Feldmann and Ronnie Kosloff, Phys. Rev. E 70, 046110 (2004). 

[13] G. Lindblad, Comm. Math. Phys. 48, 119 (1976). 

[14] Tova Feldmann and Ronnie Kosloff, Phys. Rev. E 61, 4774 (2000). 

[15] V. Gorini and A. Kossakowski, J. Math. Phys. 17, 1298 (1976). 

[16] R.S.Ingarden and A. Kossakowski, Ann. Phys. 89, 451 (1975). 



10 



[17] M. Kumakura, Y. Shirahata, Y. Takasu, Y. Takahashi and T. Yabuzaki, Phys. Rev. A 68, 
021401(R) (2003). 

[18] M. O. Scully, M. S. Zubairy, G. S. Agarwal, H. Walther , Science 299, 862 (2003). 



11 



